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Abstract 


We study the fractal structure of space-time of two-dimensional quantum gravity 
coupled to c = —2 conformal matter by means of computer simulations. We find 
that the intrinsic Hausdorff dimension du = 3.58 ± 0.04. This result supports the 
conjecture dn = —2ai/a^i, where is the gravitational dressing exponent of 
a spinless primary held of conformal weight (n J- 1, n J- 1), and it disfavours the 
alternative prediction dn = 2/|7|. On the other hand (/”) ~ for n > 1 with 
good accuracy, i.e. the the boundary length I has an anomalous dimension relative 
to the area of the surface. 
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1 Introduction 


We still do not understand the theory of quantum gravity. In four dimensions it 
has been difficult to reconcile quantum theory and gravity. We have been more 
successful in two dimensions. Matrix models 0.0.0 and Liouville held [|, ^ theory 
allow us to understand a lot about the interplay between two-dimensional quantum 
gravity and conformal held theory. 

It is sometimes said the two-dimensional quantum gravity will tell us little about 
four-dimensional quantum gravity since there are no gravitons in two dimensions. 
However, many of the conceptional problems of a theory of quantum gravity remain 
the same in two and four dimensions. In a certain way two-dimensional quantum 
gravity is as “quantum like” as a theory can be: In the path integral we perform a 
summation of all geometries with weight oneQ. This means that there is no classical 
background space-time around which we expand and our understanding of space- 
time in such a quantum world could indeed contain very important messages of use 
in higher dimensions. 

From this point of view it is somewhat annoying that precisely the structure of 
space-time is the least understood in two-dimensional quantum gravity. The recent 
introduction of the so-called transfer matrix ^ allows us to analyse in a satisfactory 
way the fractal structure of space-time for pure two-dimensional quantum gravity 
1^, 1^, 1^ and it highlights the fact that the dimension of space-time in quantum 
gravity is a dynamical quantity. Even if our underlying theory is two-dimensional 
we cannot be sure that this is the case for the quantum average. In fact, the fractal 
dimension of pure two-dimensional quantum gravity is four ! 

While the transfer matrix technique works perfectly for pure two-dimensional 
quantum gravity it is difficult to implement in the case of where matter helds are 
coupled to two-dimensional gravity. The reason is simple: when performing the sum 
over intermediate states, we should not only sum over all geometries of a certain 
kind but also over intermediate states of the matter helds. This has until now been 
impossible and another strategy has been followed: the concept of intermediate 
states is redehned in such a way that the summation can be performed [^]. The 
price paid is that the concept of intermediate state looses its direct link to the 
geometry present at the underlying manifold. To be more specihc, let us consider 
two-dimensional quantum gravity coupled to a c = 1/2 conformal held theory. At a 
constructive level, this theory is realized as the theory of dynamical triangulations 
with Ising spins placed at the centers of the triangles and the coupling of the Ising 
spins at a certain critical value. The geometry is now well dehned via the dynamical 
triangulations and we can discuss the propagation of a one-dimensional universe with 
length ii to another one-dimensional universe with length £ 2 , the two separated a 
geodesic distance D. In the case of pure gravity this problem is solved by the 
transfer matrix method which reduces the calculation to a summation over successive 
amplitudes between universes where the one-dimensional boundaries are separated 
an inhnitesimal distance. The only technical problem in the case of pure gravity 

^Except for a cosmological constant term, which has no dependence on the derivatives of the 
metric. 


2 



is that the summation over all intermediate length i of the boundaries has to be 
performed. This problem was solved in |p. When we have Ising spins at the centers 
of the triangles we have in addition to sum over possible spin assignments at the 
boundaries of length £. Presently this cannot be done analytically. This problem can 
be avoided by redefining what is meant by the distance relating two one-dimensional 
loops by only considering deformations from loops with all spin aligned to other 
loops where all spin are aligned (the prescription can be made precise). However, 
such loops might not have a well defined geodesic distance. Nevertheless, consistent 
scaling relations can be defined in terms of this modified distance. If we assume 
that the modified distance is proportional to the real geodesic distance when the 
quantum average is performed, we get the prediction that the fractal dimension of 
space-time for two-dimensional quantum gravity coupled to a conformal held theory 
of central charge c is 

= I / M , ( 1 ) 


i7(c)r 


where the string susceptibility is given by 


7(c) = 


12 


( 2 ) 


In particular, for an (m, m-|-l) conformal held theory one obtains dn = 2m, and for 
the Ising model, which corresponds to m = 3, one gets dn = Q- On the other hand, 
for c = —2 (a non-unitary conformal theory) we get 


dnic = — 2 ) = 2 , 


(3) 


i.e. the fractal dimension dn is equal to the dimension of the underlying manifold. 

An alternative prediction of dn is obtained by use of the dihusion equation in 
Liouville theory ||T9|: 


dff — 


^ Oil ^ a /25 — c -|- a /49 — c 

CX—l 


^/25^+^/T^ ■ 


(4) 


The origin of this equation is to be found in the analysis of the dihusion equation 
in Liouville theory IQ and is based on the observation that for random walks on a 
two-dimensional manifolds of area A we expect 


dim[ <^r^(f)^^ ] = dim[A^/'^^ 


(5) 


where r{t) is the geodesic distance of the random walk at the (hctitious) dihusion 
time t and the average refers to the functional integral over geometries and matter. 
In Liouville theory one can use the De-Witt short distance expansion of the heat 


kernel in terms of geodesic distance to deduce |19| that 


dim[(^r\t))j = dim[ ; 


( 6 ) 
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In eq. (|^) a_„ denotes the gravitational dressing of an (?7,+l,?7,+l) primary spinless 
conformal field, i.e. 

Jd‘^^Vg^n+i{g) j for ^^^(0 = 


where gfj,u{0 is the background metric and ^n{g) satisfies 
The requirement that (^) is a (1,1) conformal field fixes 




2n 


0^77, 


1 + •d(25-c-24n)/(25-c) 


(7) 


It is an important assumption in this derivation that it is legal to commute the 
asymptotic De-Witt expansion of the heat kernel and the functional integral over 
geometry and matter. For c = 0 one obtains dn = 4, in agreement with the 
transfer matrix prediction, while for c = 1/2 (corresponding to the critical Ising 
model coupled to gravity in the transfer matrix formulation) one obtains dn = 
(v^ -I- 7)/4 = 4.212_ For c = —2 one obtains 


dnic = — 2 ) = 


3 + VT7 


= 3.561. 


( 8 ) 


Except for pure gravity the two predictions disagree. It is the purpose of the 
present work to test if any of the two predictions is consistent with numerical simula¬ 
tions. The available analytical methods have build-in assumptions. No assumptions 
have to be made by a brute force numerical simulation of the system. Of course 
numerical methods have other problems, most notably that of accuracy. Until now 
the numerical simulations have been concentrated on systems with c > 0, mainly 
the Ising spin coupled to gravity (c = 1/2) and the three-state Potts model coupled 
to gravity (c = 4/5) 0.00 . The results have so far not been able to support 


the prediction (Q) (i.e. dn = G and 10 for the Ising and the three-state Potts model, 
respectively). However, it could be argued that these dimensions are so large that 
it would be very difficult to observe them in numerical simulations with the present 
size of lattices. It is natural to require that one should be able to probe lattice 
distances r such that 

l<r<iV^/'^«, (9) 

where N is the number of triangles or vertices in the triangulation. For systems 
with N < 10® it is clearly problematic to fulfill (^). It is however possible to 
measure the critical indices of Ising and three-state Potts models coupled to quantum 
gravity with good precision Hi- Since many of these critical indices come from 
integrated two-point functions which also measure dn it is not easy to understand 
how the critical properties come out right if (|^) is never satisfied. Moreover, it was 
recently found that correlation functions defined in terms of geodesic distance scale 
consistently with the theoretical critical indices if and only if dn ^ 4: M, O, [13 


These are indirect arguments in disfavour of (P, but of course not conclusive since 
numerical peculiarities could conspire and still allow us to determine the critical 
exponents without (P ever being fulfilled. However, in order to avoid this discussion 
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completely it is convenient to turn to the c = —2 system. This system is in a 
way the simplest coupled gravity-matter system. Within the context of dynamical 
triangulations it was the first system which could be explicitly solved 0, apart from 
pure gravity itself. Quite recently it has been possible to construct explicitly the 
two-point function using the transfer matrix methods [T^ and thereby generalising 
the results known for pure gravity |]^. 

The advantage of the choice c = —2 is two-fold. The prediction from the transfer 
matrix formulation is dn = 2, while it from Liouville diffusion is 3.562 (see (^) and 
(ID). These values for dn are so small that one has no problem satisfying (j^) for 
the sizes of systems available on the computer. Further, c = —2 is special since 
one does not have to perform Monte Carlo simulations in order to generate lattice 
conhgurations with the correct weight [|T^. As we will review below there exists a 


recursive and very fast algorithm which allows us to generate directly independent 
triangulations. In this way one can use larger systems and obtain better statistics. 
This method was first used in to investigate the fractal properties of space-time 
for the c = —2 system coupled to gravity. It was the first numerical confirmation 
of the fractal structure of quantum gravity in two dimensions. In this work the 
emphasis was put on the use of very large systems in order to get an unambiguous 
identification of continuum observables. Since then it has been shown that (a): 
finite size scaling is by far the most powerful tool for extracting critical properties 
in two-dimensional quantum gravity (i.e. the situation is similar to the one for 
ordinary statistical systems) and (b): the fractal properties of space-time have an 
interpretation as critical indices associated with two-point correlators, precisely as 
in ordinary statistical field theory |]^, |^. In this article we will show that (a) and 
(b) together with the powerful technique of recursive sampling available for c = —2 
coupled to quantum gravity makes it possible to determine dn with a precision not 
known before. 

The rest of this paper is organised as follows: In section 2 we discuss the model as 
well as the observables and their scaling. Section 3 outlines the numerical methods 
used. Section 4 contains a summary of the numerical results obtained, while section 
5 estimates the finite size shift of geodesic distance from the theoretical point of 
view. In section 6 we discuss the results and the implications. 

The results strongly support (H), i.e. it lends support to the prediction (^. How¬ 
ever, we also hnd a surprising new scaling law for the boundaries in two-dimensional 
quantum gravity. A short report of some of the results discussed in this paper has 
appeared in m- 


2 The model 

2.1 The partition function 

Within the framework of dynamical triangulations the conformal theory of c Gaus¬ 
sian fields coupled to quantum gravity is described by the following partition 
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function 

I[dxi exp-(y2{xi-Xj)^) 6 ( 10 ) 

TjveTjv '' i=i ^ \ij) i 

where TJv denotes the set of triangulations of fixed topology (which we always as¬ 
sume is spherical) constructed from N triangles. Stj^ is a symmetry factor. The c 
independent Gaussian variables xf can be viewed as placed at the center of triangle 
i. They interact with the Gaussian variables at the neighbouring triangles and J2{ij) 
denotes the sum over all such pairs of triangles. 

The Gaussian integration can be performed and one obtains (up to a constant 
of proportionality) 

Zn = Y1 T— ('^61 5 (11) 

Tjv 

where Ctj^ is the so-called adjacency matrix of the closed 0^-graph (p^iTpf) dual to 
Tat. From graph theory it is known that det^Grjv is equal to the number of rooted 
spanning trees in the graph 0^(TAr). Eq. ( pdj) serves as a definition of a model when 
c is not a positive integer, in particular when c = —2. The string susceptibility can 
be calculated in this model and it agrees with the continuum calculation in Liouville 
theory for a c = —2 theory. 

It is seen that Eq. (pdj) is special if c = —2 since in this case we can use the 
fact that det'Grjv is the number of spanning trees of 0^(T/v), i.e. the number of 
possible ways to cut the 0^(T/v) graphs of spherical topology into tree diagrams. The 
triangulations in Tpj are in one-to-one correspondence with the connected 

planar graphs with N vertices and no external legs. This can be symbolically written 
as follows: 

= E !■ (12) 

Tjv spanning trees 

in <I>Z{Tn) 


2.2 Trees and rainbows 

Let us briefly describe the combinatorics associated with the decomposition of the 
planar graphs (/)^(T/v) of spherical topology into trees and rainbow diagrams. Let 
Tn and be the number of rooted dual tree diagrams with n -|- 1 external legs 
and the number of rainbow diagrams with n lines, respectively. Especially, we have 
Ti = Rq = 1. Here, we mark one of the legs for each tree diagram and for each 
rainbow diagram in order to break the symmetry. Since any planar closed graph 
can be obtained from a spanning tree by connecting the external vertices of the tree 
by rainbow diagrams we can write ([I^ ) as 0 (See fig. 

^ Tv ^ TkTiTmR{N+2)/2 ■ (13) 

with k-\-l-\-m = N-\-2 
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Figure 1: Partition function described by tree and rainbow diagrams. 


To calculate and is easily done as follows. The tree diagrams satisfy the 
graphical Schwinger-Dyson equation shown in £g. ^a, i.e. satishes 

n—1 

Tn = '^ TkTn-k ■ (14) 

k=l 

In the same way the rainbow diagrams satisfy the graphical Schwinger-Dyson equa¬ 
tion of £g. leading to an identical equation 

n—1 

Rn = '^ RkRn-k-1 ■ (15) 

k=0 

In order to solve the Eqs. (0) and (0. we introduce the generating functions for 
the tree diagrams and the rainbow ones as 

OO OO 

Tiz) = Ri^) = E RnZ^- (16) 

n=l n=0 

Then, Eqs. ([T^) and (|^ are written by using the generating functions as 

T{z) = l + zT{zf, (T(0) = 1), (17) 

R{z) = 1 + zR{zf, {R{0) = 1). 

The solutions of (p!7D are 

= = ( 18 ) 
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Figure 2: Graphical representation of Schwinger-Dyson equations for (a) tree dia¬ 
grams, (6) rainbow diagrams. 


Therefore, one hnds 


Tn — Rn-i — 


(2n-2)! 


(19) 


n\{n — 1)! 

Using the relation (|T^, the partition function (|^ finally can be written as a simpler 
expression. 


Zn = 


N + 2' 


rFM+ii?(Ar+2)/2 


( 20 ) 


2.3 Observables 

We define the fractal structure of quantum gravity in the following way. Let us fix 
the space-time volume V. The average volume Sv{R) of a spherical shell of radius 
R is then 

Sv(R) = - -R), (21) 

where Z{V) is the partition function of gravity coupled to matter, with space-time 
constrained to have volume V, fy symbolises that the integration of metrics fulfilling 
the same constraint, denotes an arbitrary marked point and R>g(^, ^o) fhe geodesic 
distance from the marked point .^o to measured with respect to the metric g. We 
define the fractal dimension (or intrinsic Hausdorff dimension) dh of the space-time 
by 

Sv{R) ~ for i? ~ 0. (22) 

It is important to notice that the limit i? —0 is taken after the functional average 
is performed. Had we taken the limit i? —0 before the functional average we would 
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of course have obtained the result ^^dh'’= 2 since each manifold is two-dimensional. 
However, the limit i? —> 0 does not commute with the functional integral. It turns 
out that no matter how small R is there will always be numerous metrics (i.e. a set of 
metrics of non-zero measure with respect to Rig]) with the property that geodesic 
spheres of radius R consist of many connected components. For such geometries 
we cannot necessarily expect a growth as slow as R. The precise growth of Sv{R) 
for small R becomes a subtle question of entropy of different metrics and from this 
description it is obvious that the phrase “fractal dimension” is quite appropriate if 
dh > 2 . 

In the case of pure two-dimensional gravity it is a remarkable fact that one 
can calculate Sv{R) analytically 0.0 (it can be expressed in terms of certain 
generalised hypergeometric functions). One hnds 

SviR) = R^fiR/V-^), (23) 


where /(O) >0 and f{x) ~ for large x. It is seen that a dimensionless scaling 

variable R/V^^^ appears. For a general model such a dimensionless scaling variable, 
will dehne another intrinsic Hausdorff dimension d^- From (P^ and (^3]) 
we deduce that dh = dn = 4: in the case of pure gravity. For general model we can 
write (p^) as 

Sy{R) = (24) 

where 

Fi{x) ~ X -C 1, (25) 

and Fi{x) goes to zero as 


for X going to inhnity. Eq. (|2^) has the form 
of a typical finite size scaling relation and we expect it to be valid not only for pure 
gravity, but also for gravity coupled to matter. Since Sv{R) is easily measured in 
numerical simulations we can use Eq. 

(El) and (E5|), one finds 


and (|25|) to extract dh and du- Using 


Sv(R) ~ -e. 1. 


(26) 


If space-time for large V has the same fractal properties at all scales, one expects 


dh — du- 


(27) 


However, in our numerical simulations we do not assume this property dh = dn- It 
is one of our purposes to check if eq. (pTl) is realized for the c = —2 model. 

Let us briefly describe how the above continuum description translates to the 
framework of dynamical triangulations. To a triangulation T/v we can unambigu¬ 
ously associate a piecewise linear manifold with a metric dictated by the length 
assignment £ to each link. From a practical point of view we use instead a graph- 
theoretical distance between vertices, links or triangles. In the limit of very large 
triangulations we expect that the different distances when used in ensemble averages 
will be proportional to each other. To be specific we will in the following operate 
with a “link distance” and a “triangle distance”. The link distance between two ver¬ 
tices is defined as the shortest link-path between the two vertices, while the triangle 
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distance between two triangles is defined as the shortest path along neighbonring 
triangles between the two triangles. In this way the triangle distance becomes the 
link distance in the dual (f)^ graph. 

In the following we will report on the measurement of quantities related to the 
fractal structure of quantum space-time: The total length (/) and the higher mo¬ 
ments (/"■) of spherical shells of (geodesic) radius r, and the distribution function 
p(r, 1) which measures the (average) number of connected components of the shells 
of length I and radius r. 

More precisely, let us consider the class of triangulations Tjy which are dual to the 
connected closed 0^-graphs. The number of triangles (or vertices in the 0^-graphs), 
N, plays the role of volume. If £ denotes the link length of the triangles the relation 
to the continuum volume is Id ~ e'^N, and we want to take a limit where V is fixed 
while £ —> 0 and N to infinity. We consider a spherical ball of radius r and its 
shell for a given triangulation T/v. The spherical ball consists of all vertices with 
link distance r' <r and the spherical shell consists of all vertices with link distance 
r, where the distance is measured from a given vertex vq which is considered as 
the center of the spherical ball. In the same way we can define the spherical shell 
in terms of triangle distance. We will use both definitions in the following and we 
expect that after taking the statistical average they will be proportional to each 
other and that they will not affect the universal properties of correlation functions 
. The spherical shell in general consists of a number of connected components if 


we define a connected component of the shell of vertices as a maximal set of vertices 
in the shell where all vertices can be connected via links in the shell. If we take 
the average over all positions of vq and all triangulations T^r, we get a distribution 
P]\f{l,r) of the length I (measured in link units) of the connected components of the 
spherical shells of radius r, i.e. 


CO 


1=1 


(28) 


In particular we introduce the special notation njsiir) = and since nj^ir) is 

the discretized version of Sv{R) we expect the fractal dimension to be related to 
nAr(r) by 

n^(r) ~ 1 < r < . (29) 


According to the general scaling arguments mentioned above 
the following behaviour for nN{r): 


[0, m o 


we expect 


nN{r) ~ 


r 

]\[l/dH ’ 


(30) 


and we expect Fi[x) to behave as ^ for small x and to fall off rapidly when 
a; > 1. 


3 Numerical method 
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3.1 The recursive algorithm 

The recursive algorithm takes advantage of the factorisation property of the partition 
function (see Eq. (|T^) in order to construct a typical configuration of c = —2 gravity. 
One constructs a rooted 0^ tree with the correct probability and then connects the 
outer links of the tree with a rainbow diagram also constructed with the correct 
probability. According to Eq. ([l^) the branching probability to divide a rooted tree 
diagram with n + 1 external legs into two different rooted tree diagrams with k + 1 
and n — k + 1 (“1” counts the root) external legs is given by 


w(n, k) 


TkTn-k 


(31) 


If we want to construct a surface with N triangles {N must be even) we need to 
construct a tree with + 2 external legs. In practice, we start from the root, which 
has a tree with + 2 external legs attached to it. Then we proceed with branching 
the root into two trees with k + 1 and N + 2 — k external legs (remember that we 
also count the root leg), where k is computed from Eq. (|^) . At each step we assign 
the number of external legs of the tree attached to each link according to the same 
formula and we keep an ordered list of the external legs of the whole tree. We add 
one such link to the list whenever k = 1. 

Then we proceed to connect the external legs with a rainbow diagram. This is 
possible since the total probability is the product of the probability of constructing 
the rooted tree Eq. (^Tj) and the probability of constructing the corresponding rain¬ 
bow diagram. According to Eq. (|I5D , the probability of splitting a rainbow diagram 
with n lines into two rainbows with k and n — 1 — k lines is given by: 

u{n, k) = ^ ^ = w{n -|- 1, /c -1- 1), (32) 

Rn 

since Rn-i = Tn- In practice we start from the root leg of the tree diagram and we 
split the rainbow containing N/2 + 1 lines into two parts containing k and N/2 — k 
lines. Then we can connect the root leg with the appropriate member of the list 
containing the external legs of the tree mentioned above and proceed accordingly 
until we connect all external legs. 


3.2 The simulations 


The simulations are performed by generating a number of statistically indepen¬ 
dent conhgurations using the algorithm mentioned above. We use the high quality 
random number generator RANLUX ^ whose excellent statistical properties 
are due to its close relation to the Kolmogorov K-system originally proposed by 
Savvidy et.al. |^, 0 in 1986. We centered our effort for good statistics on system 
sizes ranging from 2000-256000 triangles. The number of configurations obtained 
depends on the lattice size and on the observable that we measure. We choose 
20 random vertices/triangles on each configuration in order to perform correlation 
function measurements. We need to collect more statistics to test Eq. (^), where 
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we have between 4.2 x 10® and 1.6 x 10® configurations. For the 128iF and 256iF 
lattices we have 6 x 10® and 2 x 10® configurations respectively. It was possible 
to extract useful information about the short distance behaviour of the two point 
function njv(r’) by generating a smaller number of configurations for system sizes 
having 512000-8192000 triangles. We got 10000 configurations for N = 512000 and 
1024000, 4000 for N = 2048000 and 386 for N = 8192000 and we measured nN{r) by 
choosing many more initial random triangles/vertices. It is not possible to extract 
dn by using finite size scaling with so low statistics. 

In order to measure the moments and their scaling properties we need a 

factor of 10^ less configurations: We have approximately 50000 configurations for 
each lattice size. Unfortunately, the computer effort for making the measurements 
is comparable to the one needed to test Eq. (|5I1| ) with enough accuracy. 

One subtle point in the simulations is the computation of the branching numbers 
given by Eqs. (^) and (^). Given the number n we need to compute k. We do 
this by choosing a random number r in the interval [0,1). Then, e.g. for the case of 
trees, we compute 

k 

W{n,k) = '^w{n,i) ^ (33) 

i=l 

and we choose k to be the integer such that lU(n, k) > r. Using the symmetry of 
w{n,k) around n/2 we can substantially reduce the computation time. Moreover 
w{n, k) is best computed from the recursive formula 


w{n, k) 


{2k — 3)(n — /c + 1) 
k{2{n — fc) — 1) 


tc(n, /c — 1), 


1 ) 


n 

2{2n - 3) ’ 


(34) 


and extra care must be put in the computer program so that we obtain the correct 
probabilities due to overflow. For this reason we have tested the distributions of k 
for given n obtained from the computer versus the theoretical value Eq. (^). The 
results for n = 16000 are shown in fig. ^(a) where we plot w{n,k) from Eq. (B 
together with the results obtained from using our program 10® times. A measure of 
the agreement is shown in fig. §(b) where we plot the relative deviation Aw{n, k) = 
{w{n, k) theoretical “ w{n, k) measured) /w{n, k) theoretical- In the latter case, the data is 
smoothed using Savitzky-Golay filters to interpolate the nearest 100 points to a 4th 
order polynomial. As is seen, the distribution is symmetrical and the relative errors 
average out to 0 very nicely indicating that the deviation is pure statistical noise. 
We have checked that \Aw{n, k)\ ~ Admeasurements- 


4 Numerical results 

4.1 The fractal dimension 

We have measured the fractal dimensions dn and dh in a number of ways to be 
described in the following. Since we are considering discretized, finite systems as 
approximations to continuum systems (although we try of course to stay as close to 
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continuum physics as possible) it is important to use independent ways to approxi¬ 
mate the same continuum physical observable. This will often give a more reliable 
idea of how close we are to the genuine continuum quantity since it supplements 
the systematic large N study of a single choice of discretization of the continuum 
observable, where systematic cancellation of fluctuations might sometimes underes¬ 
timate the real discrepancy between the measured quantity and the unknown value 
of the continuum observable. 


4.1.1 Short distance behaviour 


One can try to use directly the short distance behaviour 
dh- This was the method used in the pioneering work |T6 


of nN{r) to extract 
The problematic aspect 
of the method is to what extend one has to take seriously the natural requirement 
1 -C r We know now from the exact solution of pure two-dimensional 

quantum gravity that both limits have to be respected quite seriously, but that a 
so-called shift r ^ r + a helps to almost remove the requirement for the lower limit 
n]. Later we will discuss various theoretical and “phenomenological” motivations 
for this shift. Presently, let us just use it as an additional £t parameter. In table ^ 
and table 0 we have shown the results of a fit of the form 


Un 


(r) = C {r + aY^ 


(35) 


for various cuts of Umin < r < Umax- We see here a discrepancy between dh obtained 
from the link-distance measurements and the triangle-distance measurements, re¬ 
spectively. If there is a discrepancy between the measured dh for triangle and link 
distance for hnite N one would expect it to be “maximal” when we concentrate on 
small r, like here. For small r the triangle distance is quite rigid in the sense that 
each triangle has at most three neighbours, while a vertex can have any number of 
neighbouring vertices. This rigidity is also reflected in the fact that we have to use a 
much larger shift a for triangle-distances. There is a weak tendency for the link-d/i 
to decrease with N while the triangle-d/j shows a similar weak tendency to increase 
with N. From this one would conservatively estimate, assuming that they have a 
common large N limit, that 

3.53 <dh< 3.60. (36) 

In hgs. ^ and ^3 we have shown the behaviour of short distance behaviour of nhj{r) 
for the two distance measures and the whole range of N. 

We will now show that the discrepancy between Imk-dh and triangle-d/j indeed 
decreases when we use the hnite size scaling (|30D to invoke the whole range of r and 
also that we can get a much better determination of dh- 


4.1.2 Collapse of distributions 


Before using the hnite size scaling relation ([30| ) let ns motivate the use of the “shift” 
a which is essential for obtaining high precision results. The need of this parameter 
is well known from earlier studies of conformal held theories with c > 0 coupled 
to quantum gravity |11, 13|. One obvious, “phenomenological”, motivation for this 
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shift parameter (which we have just copied from standard hnite size scaling theory) is 
as follows: for hnite N we expect some discrepancy compared to continuum results, 
typically parameterised by “the number of points” L corresponding to the linear 
size of the system. In particular we can write 

R r const. const. 

^ “ yi/dH ^ jyi/dH I ^ L2 

or, since is precisely a typical measure for the linear extension of the system: 



R 


r + a const. 
jyl/dH (^jyl/dny • 


(38) 


The parameter a, which is considered a shift in r, incorporates the hrst order correc¬ 
tion. Below we derive a theoretical value for a in the case of pure gravity and where 
we use triangle distance. However, here we consider it as a purely phenomenological 
parameter, which in principle can be different for different scaling variables and will 
be different if we use different distance measures (viz. link and triangle distances). 

The raw measurements of n]^{r) produces distributions with N ranging from IK 
to 512i^. We can now try to ht them to the scaling ansatz (^Ol). We have two 
parameters available, dn and a. 

The collapse of the distributions is performed by using two different methods. 
The hrst one is identical to the one proposed in and used with great success 
also in BUI- One makes a non linear ht of the form Pn{x)e where Pn{x) is a 
polynomial of order n in x, to the rescaled according to Eq. (|30D distributions nj^{r) 
for a given set of lattice sizes {N}. The dn) per degree of freedom is computed 
for a given set of parameters a and dn- For each value of the shift a, the optimal 
value of dn is computed from the position of Xmin(®)- 

We also used a second method which has the advantage of being much faster. 
We should also note that it is also quite successful even for very small lattices 
{N < 500) where the hts used in the hrst method fail to yield reasonable results. 
For a given set of lattice sizes {N} and parameters (a, du) we compute a cubic spline 
interpolation to the rescaled distribution for each value of N . x^(0) dn) is computed 
by adding in quadrature the distances of each point of the other distributions from 
the interpolation function reweighted by their errors and properly normalised to 
correspond to a per degree of freedom. The best values for dn computed this 
way are identical to the ones computed using the hrst method, although x^(0) du) is 
slightly smaller and steeper yielding smaller errorbars (~ 10-30%) for the computed 
quantities. In this paper we report the larger errors computed from the hrst method. 

In table ^ we have shown the results of the best hts both for the link-distance 
distribution and the triangle-distance distribution. The hts are divided into groups: 
the upper one determine du and a by comparing distributions nN{r) for successive 
pairs of A^’s (16i^ and 32i^, 32i^ and 64i^, etc.) In this way the [dn-i a) dependence 
on N becomes clear. The middle and lower groups joins three, respectively four 
successive N^s in the determination of dn and a. In particular for the top group 
we see that du has a clear dependence on N\ If we use link distances dn^N) is 
systematically decreasing, while dniN) is systematically increasing if we use triangle 
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distances. Under the assnmption that the trend will continue for larger systems and 
that the two distance measures are proportional to the genuine continuum distance 
in the scaling limit we conclude that 


3.55 <dH < 3.61. 


(39) 


These values of dn (and the appropriate values of a) actually yield excellent hnite 
size scaling for the whole range of N is illustrated in fig. ^a where we have shown 
fW 

{x) for the whole range of N from 2K — 256K. The importance of the inclusion 
of the shift a if we use the triangle distance is illustrated in fig. |b where we have 
plotted the functions obtained from n]s[{r) with the same choice oi dn but 

with a = 0. The same plots are shown in fig. 0a# in the case where we use 
the link distance. Here a is smaller and not of the same visual importance as for 
triangle distances. However, in the actual hts they play an important role for the 
link distances as well if we want to extract consistent values of dn with acceptable 
values for the hts. In £g. || we see that the value of dn depends strongly on a. 
However, we can use the systematic behaviour of dn as a function of N and the 
shift parameter a to obtain a better estimate of dn than the one provided by 
In £g. ^ we have shown the complete Xmin(®) which was used in the top group of 
table 1^. Each point on the two subhgures represents a specihc choice of the shift a. 
For a given N we now End the value dnia, N) which minimises the of fhe collapse 
of the njsiir) and nj^Fir). In table ^ we recorded just the minimum as a function 
of a and the change of the minimum as a function of N is seen quite clearly in £g. 
1^. However, let us for each N plot the (a, N) of £g. as a function of the shift 
a. This is shown in Eg. H. We get a number of straight lines which with very good 
accuracy intersect for one value of oq. The simplest phenomenological explanation 
of this fact is that Oq is the correct value from (^), since this implies that 


dn^a, N) = duiao, oo) + 


ddniao, oo) 
da 


+ 


d'^dniao, oo) 1 
daod{l/N) N 


(a 


Oo) + ■ ■ ■. 


(40) 


From this argument the correct inEnite volume limit of dn is determined from the 
Egures up to the accuracy with which the curves actually cross in a single point. 
This interpretation conErmed by the fact that both the triangle and link dH{a,N) 
curves cross at the same value of dn, but for very difierent Oq. 

We conclude from the data (see figure) that 


dn = 3.57 ±0.01. 


(41) 


4.1.3 Average radius 

In this section we will use the average radius of a universe to extract the intrinsic 
Hausdorfi or fractal dimension. From the definition of ni^{r) the average radius of 
universes with volume N is 

1 CO 

{'f') N = ^N{r) ^ . (42) 

r=0 
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Obviously, ( ^^j ) could itself serve as a natural definition of dn- By measuring n^ir) 
we can record (r)^ as a function of N and hence determine dn- As usual we have 
to introduce the shift a in order to account for lowest order discretization effects. 
Now, let us dehne 

(r+ a)^ 


Ra,N{d) — 




(43) 


We determine the values of a and dn in the following way: hrst we measure (r)^. 
for a certain number of different volumes W of the universes, ranging from N = 2K 
to N = 256K. For a given a we choose, for each couple Ni, Nj of iV’s, the such 
that 

Ra,Ni{d^) = Ra,Nj{d^) . (44) 

For this choice of W, Nj we bin the data and estimate an error Then we 

determine the average 


dn = 




and compute 


X^(a) = 


# pairs ^ 

1 (d*^ - dnf 


T. 


(45) 


(46) 


# pairs ^ (Sd'iy- 

The preferred pair (a, dnia)) is determined by the minimum of x^(a). This method 
works quite impressively. In £g. [lO| a and fig. |lO| b we have shown the intersection of 
the curves i?a,Af(d) as a function of d for the optimal choice of a for link-distance 
and triangle-distance measurements, respectively. The important points are that in 
both cases there exists a value of a where all the curves intersect with high precision, 
and that the range of a where x^(®) is acceptably small, i.e. 0(1), is quite small. 
Hence d^ is determined with high precision. In Fig. |T^ and |^b we show x^(®) 
link-distance and triangle distance measurements, respectively. In this way we get 


duiam) = 3.573 ±0.005 


(47) 


both from the link-distance measurements and the triangle-distance measurements. 
The values of Om are 


= 0.13 ±0.01, 


and 


= 5.00 ±0.05, 


(48) 


for the link-distance and triangle-distance, respectively. In table ^ we have listed 
the determination of (d/f(am),am) for various cuts in the lower values of N. The 
constancy and consistency of the results are truly remarkable. 

In (^71) we have estimated the error as follows. Define an interval of acceptance 
[amin, Omax] of o by demanding that x^(®) < ‘^X* where x* = niax{l, x^(am)} and 
hnd the variation of d(a) in this interval. After this we repeat the whole procedure 


by making various cuts in the pairs of W’s included in (^Sj) and (^), discarding 
successively the smallest W’s. 


We note that the values (|47|) are in perfect agreement with (^T]) . We have now 
four independent measurements (two different methods, two different dehnitions of 
length) which all give dn = 3.57 ±0.01. 
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4.2 Boundaries 


We now turn to the measurements of These observables are constructed from 

PAr(/,r), which can readily be measured in the simulations. We recall the situation 
in pure two-dimensional quantum gravity which can be solved explicitly and where 
dh = 4. In this case we have (for small r): 


n>l. 

(49) 

4 in this case, i.e. we obtain 


(niTN-i^y, n>l, 

(50) 


as one naively would have expected for a smooth two-dimensional world. Only the 
hrst moment behaves anomalous, as it has to do if dh 7^ 2. From these c = 0 
considerations it is unclear what to expect for c = —2 for the higher moments. If 
dim[iV] = dim[/^], then from scaling arguments, we expect 

HXn ~ , X = ^ . (51) 

However, our measurements are consistent with the following scaling relations 

{nr^N ~ for n > 2, (52) 

which implies that dim)/”] = dim[r^"’] for n > 1. Eq. (|5^ indicates that we have 

~ for 1 < r < n>2. (53) 

We have shown this relation in Eg. |T^ for n = 2, 3 and 4 for the case where we 
use link distances. It is remarkably well satisfied. In table ^ we have shown the 
more detailed result of this short distance analysis. It should be compared to the 
short distance analysis presented in table (and £g. Pd) for the hrst moment. The 
behaviour of the exponents 2n for the higher moments as functions of N are clearly 
more consistent and systematic compared with the behaviour of dniN) for the hrst 
moment. In addition one can use smaller values of r and there is not the same crucial 
dependence on the shift a as for the hrst moment. All this points consistently to 
smaller short-distance discretization ehects for the higher moments. Note also that 
the value of the shift a is diherent. 

Let us now turn to the hnite size scaling analysis of (|52D. Again we introduce 
as a hrst phenomenological correction the shift r —»• r -|- a as in (^3]) to hnd the 
best scaling function F„(x) for a suitable range of W’s. We have shown F„(x) for 
n = 2, 3 and 4 for the values of a which provide the best scaling function in hg. |^. 
The more detailed analysis is given in table Note that the results are perfectly 
consistent with the same analysis for the hrst moment (table but the shift a is 
diherent and in fact not very well determined (see hg. |T^). In principle this is a 
good thing, but it implies that we cannot use the shift in the same constructive way 
as for the hrst moment and get a high precision measurement of dn- 
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The scaling ansatz (pl|) seems to be ruled out. In fig. |^a we have shown the 
best overlap functions ^ 2 ( 2 ^) for the second moment {l‘^)rN using (|^ . It should be 
compared to fig. 


13, where we have used the ansatz 


Clearly 


IS superior 

to (ED- One may, however, consider the possibility of a dimensional relation of the 
form: 

dim[/”] = dim[r^”*'^“^^ ] . (54) 


In this case, for given dn, one can use the relation = N^n(l-e)/dH in 

order to determine the value of e. We get an upper bound on e by using the lowest 
value of for du (obtained by some of the other methods discussed above) and then 
fitting to the data. In this way we obtain e < 0.03. We consider the existence of 
such a small e 7 ^ 0 unnatural. 


4.3 Distribution function 

Let us finally turn to the measurement of the distribution function pAr(r,/), which 
provides the complete information about the moments ^ ‘^use of pure 

gravity one can calculate pAr(r, /) in the limit N ^ 00 and one finds P: 

p!Sr°\r,l) ~ (55) 

where the function (z) is 

= + ^ + ( 56 ) 

Note that this form of p^^{r,l) explains (^O]). If e denotes the cut-off (the lattice 
spacing in the triangulation) we can write 


in 


c=0 

r,N=oo 


dirp<£(r,[) 


= r^”/ d 2 2 ”G<">(j). 

Je/r^ 


(57) 


For n > 1 the lower limit has no implication for the integral and can be dropped. 
However, for n = 1 the leading contribution in the limit £ —> 0 comes precisely from 
the lower integration limit and we obtain: 


(or^=oo 





(58) 


The cut-off dependence ensures that the real dimension of (/) is equal to that of 

Vn. 

Since we have verified with good accuracy that the higher moments for c = —2 
and c = 0 has the same r dependence for large N (see (^^ and (|53[)), we know that 

PAr(r,/) ~ ^G(//r^) for N ^ oo. (59) 


The function G{z) need not be identical to G*^°^( 2 ;). In fact it cannot be identical 
since dh{c = — 2 ) ^ dh{c = 0). Since we expect that the origin of the different 
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behaviour of the hrst and the higher moments is the same for c = 0 and c = —2 we 
know that the small ^ behaviour of G{z) must be 


G{z) ~ 


1 


dh ^ 3.57. 


(60) 


Fig. |l^ a displays that the link-distance distribution function r‘^p^{r,l) is to a 
good approximation only a function oi z = l/r‘^ as long as < 1. For small 

values of 2 ; the log(G(^))-curve is approximately a linear function of log( 2 ;). From 
(|60|) we expect 


logG(z) ~ 


+ 1 


\ogz 


for log 2 : —> —cx). 


(61) 


A measurement of the slope gives dh = 3.60 ± 0.03 in good agreement with short 
distance behaviour (Or given in table for the link distance. We emphasise 
that this agreement is no surprise. There has to be agreement. However, as we have 
seen, some treatments of the data set can produce very precise measurements of dh. 
The slope of logG(z) does not belong to this class. 

Finally, in £g. IHc, we have shown r‘^pN{f', 0 ^ somewhat larger z interval. We 

know from measurements of the large r behaviour of the moments that they 

fall off fast when x = > 1 (see (^) and fig. |l^-0c). The same behaviour 

has to be coded in pN{r,l). The simplest guess: 


PN(r,l) ~ ^G(z)H{x) 


(62) 


is not satished since this would imply that the scaling functions Fn{x) were pro¬ 
portional to x‘^^H(x), which is not very well satisfied numerically. However, the 
ansatz ([6^ ) is not too far from the truth either. Since the curves shown in £g. ^c 
are very similar to the ones for pure gravity one can test that an ansatz like (|6^, 
replacing G{z) with and H{x) with the function f(x) from (P3|). The ansatz 

will produce curves qualitatively very similar to the ones shown in £g. ^c. Indeed, 
the curves consists of three major parts. Let y = log//r^. The right part of the 
curve goes as —and comes from the exponential decay of G^^\l/r‘^) for positive 
y. The “straight-line part” is —{5/2)y which dominates for negative y, as discussed 
above. The “bump” is the joining of these two asymptotic regions. Finally, for hxed 
llj\j2/dH large negative y we hit the region where f{x) goes to zero. It will 
contribute with a term (for dn = 4) 



h 


(63) 


This term depends on ljN‘^1'^^ and consequently the decay of p^ir,!) for large neg¬ 
ative y will depend on 1 /N'^Fh precisely as is shown for the c = —2 case in £g. 


16c. 
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5 Theoretical estimate on the finite size shift of 
geodesic distance 

In eq. (|5^ ) we introduced the shift as part of a general hnite size expansion. Since 
it plays an important role in the hts it is worthwhile to provide a theoretical under¬ 
standing of the magnitude of the shift. We shall limit ourselves to an analysis of the 
c = 0 model where it is possible to perform a theoretical analysis at a discretized 
level. There are other technical differences between the theoretical and the numerical 
analyses which will be explained in order. Thus the comparison will be qualitative. 
Yet we believe that the present analysis provides a theoretical conhrmation of the 
existence of a hnite size shift of the geodesic distance. 

The hnite size shift a will be shown to incorporate the ehect of next higher-order 
correction in a lattice spacing parameter. Here, we denote the two-point functions 
at the discrete level and at the continuous level as Gg{r) and Gg{R), respectively. 
The two-point functions are related with Sv{R) and njsiir) as 

Gg{R) = [ dVe->^^VZ{V)Sv{R), Ggir) = NZ^n^ir), ( 64 ) 

N 


where we have introduced g a.s a coupling constant of matrix model, r as the geodesic 

distance in the discrete level, /i as the cosmological constant, and R as the geodesic 

2 

distance at the continuous level, g and r are related with /r and Rhj g = gce~°“^ ^ 
and R = f3re^'^ ^ where a and (3 are constant, and e is the lattice spacing parameter 
in the triangulation. 

According to the arguments in ref. 0], the continuum limit of the two-point 
function is 

Gg{r) = e‘^^^-^>{G^{R) + e^^G^^\R) + 0(8^^} , ( 65 ) 

where u and rj are constant. Suppose that the following relation 

Gp{R)^l3aZc:^{R), ( 66 ) 

is satisfied, we obtain 

Gg{r) = e^^^-^>{Gg{R + , ( 67 ) 


where this a will be identihed as the hnite size shift of geodesic distance. Thus, one 
can incorporate next higher-order correction by redehning the geodesic distance at 
the continuous level by instead of by simply taking R = jdre'^'^. 

Now, let us carry out a concrete calculation. Here, we restrict to the analysis of 
pure gravity (c = 0 model) because this is the only case where the two-point function 
is theoretically known as a function of the geodesic distance R. In dehning geodesic 
distance at the discrete level, there are two types of decomposition of triangles; 
slicing decomposition and peeling decomposition [Q. Gg{R) is known for both 
cases while the discrete two-point function Gg{r) is known only for the case of the 
peeling decomposition. In the peeling process, a triangle can be peeled off with 1/1 
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step forward of geodesic distance (Z: loop length at the discrete level) while fnll one 
step forward of geodesic distance is realized in the one slicing process. 

In the peeling decomposition! ^|, the two-point fnnction is related with the gen¬ 
erating fnnction of the disk amplitnde at the discrete level 0 ], 


d 

Gg{r) = —Fg{x{x,r)) 


x=0 


g dr 


Fg(x(0,r)), 


( 68 ) 


where we start from the bonndary of disk with length Zi = 1 at r = 0. Here, 
(x) = YxiX^Fgil) is the generating fnnction of the disk amplitnde Fg{l) with one 
marked link and bonndary length I at the discrete level. The fnnction x{x, r) satishes 
x{x, r = 0) = X and 

dx fg{x) dx 

dx 


fg{F) ’ 


dr 


= 9fg{x) 


(69) 


In the case of one-matrix model [^], Fg{x) and fg{x) have the following forms, 

+ fgix) = il-C2x)^J{l-Cix){l-Cox), (70) 




9 


where cq < 0 < ci < C 2 . The solntion of 
... 1 ' 


IS 


x(x,rj = 


1 - 


<5i 


‘^2 I smh'^(6or+ smh 


(71) 


where 


X _9 d W ^ _ (C2 - Ci)(c2 - Co) ^ _ Co(c2 - Cl) 

5o-^V(c2-Ci)(c2-Co), c2(Ci-Co) ’ ~ ~ rJr. - r.Y 


C2(Ci - Co)' 


Next, let us consider to take the continuum limit of (| 68 |) . The continuum limit of 
the disk amplitude is taken by imposing x = Xce~^^ and g = gc€~°"^ 

Fg(x) = (const.+ fg(C) + , 


mo = {<:-^)FFF. 


(73) 

(74) 


where the constants, bi and 62 , and the critical values, Xc and gc, are regularization 
dependent. The continuum limit of x(0,r) is 


£(0,r) = Xcjl - eCo{R + (3aiy/e) + 0(e^)^ , (75) 

where 

Qo(R) = -y/JI+^y/JIcoth^(^^fj.F^R^ , (76) 

and R = (3ry/e, which means that i/ = 1/4 in pure gravity. Substituting ([731) and 
(|75|) into (| 68 |) , we hnd 

F'gi '^) = (const.)£^'^^^^|— 62 Co(.R + Po.l^/e) + V^fgiCoiR)) + 0(^)|, (77) 
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which tells us that rj = A. Here, we obtain the following nontrivial relation from 

( 0 ). 

= -2/„(C„(fl)). (78) 

Then we hnd the property (|66D , i.e., we hnally find 

Gg{r) = -(const.)£^/^^|Co(^ + |3al^/e + + 0(e)|. (79) 

Thus, we obtain 

^modified ^ ^y2v ^ ^ ^ 2^) ' 

In the one-matrix model which describes pure gravity, the critical values are Xc = 
(3^/^ — 3“^/^)/2 and Qc = 1/(2 ■ 3^/^). Then, we find 62 = ■\/3(v^ -|- 1)“^'^^, (3 = 
(a/3 -|- 1)^/^/(2a/3), and ai = a/3. Therefore, the total finite size shift becomes 

a = Oi -|- 02 = 2a/ 3 -|- 1 = 4.46... (oi = \/3, 02 = a/s -|- 1). (81) 


We now try to evaluate the two-point function where the starting boundary 
length is 3 at r = 0. In this case we have to take the third derivative of the 
generating function of the disk amplitude instead of the first derivative of (| 68 [) , 






1 3 

+ 




x=0 

1 


(82) 


d 1 
+ 


g\fgJ dr g^\drj \ g dr 


d 


1 d 


FAi) 


x=0 


The continuum limit of 

G!;-=‘l(r) = (const.) (l+/ 3 a 3 v^^ + 0(£))Gj;—l(r), 
where G^g^^^\r) is Gg{r) of ([791) and 

m/fg)' 


03 = lim 


°£/(1//9 


x=0 


(83) 


(84) 


In the one-matrix model, the value of 03 is 03 = 3/2. For the case with the boundary 
disk length Zi = 3 at r = 0, we obtain the total finite size shift. 


Qi — Oi -|- Oj2 T 03 — 2'\/3 4” 2 — 5.96.... 


(85) 


The similar calculation can be done for the triangulation without two-folded 
links (two links are on top of each other) which come from the kinetic term in the 
matrix model. In this case the disk amplitude is 


X 


Fg{x) = ^ - 2 “ ^ + ^fgix ), fg{x) = (1 - C 2 a;)Vl-cia;, 


( 86 ) 


22 











where ci = A/i^gc^) and c\ — C 2 Ig^ = —9>lg. Since fg{x) in (PB| ) is formally the same 
as fg{x) in (|70D with the snbstitntion cq = 0 , x{x, r) of the present model is formally 


the same as x{x,r) in ([7l|) with cq = 0. In this case, the continunm limit is realized 
around the critical values, Xc = 1/(2 ■ and gc = 1/(2 ■ 3 ^/^). The continuum 
limit of Fg{x) is ([ 73 |) with 62 = 2 / a / 3 , and that of a;(0,r) is (^) with (3 = l/(2v^) 
and Oi = 2 ^/ 3 . Then, for the present model we find 

3 3 

a = Oi + 02 = 2'/3 + 2 ~ 4.96... (oi = 2 \/ 3 , 02 = -). 


(87) 


In the present case with the boundary disk length /i = 3 at r = 0, the total finite 
size shift is 

16^3 3 „ , 

o — Oi + 02 + 03 — —p-1" Q ~ 7.04..., ( 88 ) 

0 Zi 

where 03 = 6v^/5 which is derived from (| 8 ^) . 

In the slicing decomposition in pure gravity P, we have 02 = \/(38 + 21v^)/6, 
because of the property that (|78[) is independent of regularization. However, we 
failed to obtain oi because we do not know the concrete expression of x{x,r) in the 
slicing decomposition. 

The above arguments show that the hnite size shift o depends on the detailed 
prescription for the starting point and the class of triangulations used. As a further 
source of ambiguity we can mention that the counting of distance in the theoretical 
considerations above and in numerical simulations differ in the following way: In 
our numerical simulations we first mark one of the triangles on the sphere, and then 
measure the distance from the marked triangle. The marked triangle might be a 
normal triangle or a triangle which corresponds to a tadpole in the dual lattice. Thus 
the marked point is located at the center of the marked triangle in our simulation 
while the boundary of the marked triangle is the zero geodesic region in the above 
theoretical analysis. In £g. we describe the definitions of the distance both in the 
above theoretical analysis and in our numerical analyses. Therefore, the geodesic 
distance in the numerical simulation r®™ is related with the geodesic distance r as 
7 -sim ^ x/ 2 . Since _R™°dified _ _ p(j.siTa _|_ ^sira^piu^ which reproduces 

the universal function we hnd a®™ ~ a — 1 / 2 . Another slight difference 

between the theoretical and numerical determination of the shift a is the different 
definition used for the boundary lengths. In the numerical simulations we count the 
number of triangles at a given distance r®™ whereas in the theoretical approach we 
count the number of links at distance r. Finally, in the computer simulations we 
use a modified version of the slicing decomposition. This implies that we cannot 
directly compare with the theoretical results derived above which use the peeling 
decomposition to estimate the geodesic distance r. However, when we compare with 
the numerical simulations with c = 0 it is found that the shift is4<a<6|ll 

and (^. 


which is consistent with eqs. ( 0 ),( 

For more general models where c 7 ^ 0, we have not succeeded in showing the 
property (| 66 |) which is used in proving the existence of the shift a. On the other 
hand, the numerical simulation in c 7 ^ 0 models strongly supports the existence of 
the shift a of the same order of magnitude as for the c = 0 model. 
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Figure 3: (a) The distance used in the theoretical analysis, (b) The distance used 
in our numerical simulation. The difference of the distances between (a) and {b) is 
about 1 / 2 . 


Finally, we consider the shift in {l^)r- In the peeling decomposition, {l"')r is 
described by 

_ d r d 


in 


dx V dx 


Fg{x{x,r)) 


x=0 


(89) 


where we start from the boundary of disk with length /i = 1 at r = 0. Note that 
from ([75|) we obtain 

d \ d 

x^ = --^ + 0{e). (90) 

edC,o 


dx 


Substituting (p0[) into (| 8 ^ and carrying out the similar calculation as before, we 
find for n > 2 , 


(r), = (constOe’ ”^(^) {/(.(Co(fi+/3«iv^)) + («)nst.)y?(Co(fi))^+0(i;)}' 

(91) 

For n > 3 the finite size shift comes only from Oi because 02 is absent. When we 
start from the boundary of disk with length /i = 3 at r = 0, the shift is 01 + 03 . 
The values of o for {n > 3) are all equal but are smaller than that of o for 
{l)r by 02 . On the other hand, for n = 2 the property (| 66 |) is broken, because Cq 
is not proportional to {d/dR)fg{(o). So, we cannot expect a clear shifting property 
for {R)r in the numerical simulations. 


6 Discussion 

In this article we have taken advantage of the special recursive sampling possible for 
the d = —2 theory coupled to quantum gravity. The quality of the numerical data 
obtained this way is much better than the quality of data obtained with a similar 
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computer effort by ordinary Monte Carlo simnlations. The pnrpose of the present 
work has been to measnre the fractal strnctnre of space-time in the c = —2 theory 
conpled to qnantnm gravity with a precision which has not been available before, 
by combining the high qnality data with the techniqne of hnite size scaling. 

In this way we obtained, with a conservative error estimate, 

dn = 3.58 ± 0.04, dh = 3.56 ± 0.04, (92) 

in perfect agreement with the theoretical prediction 3.561 ■ ■ ■ from the diffnsion in 
Lionville theory (see (^) and (|^)), and in disagreement with the prediction given by 

®-(®. 

Thns eq. (§) is strongly favonred as the correct formnla for the fractal dimension 
of space-time in the case where the matter helds conpled to gravity have c < 0. 
In this region eq. (^ has nicer properties than eq. O. since one naively expects 
that dn^c) ^ 2 for c —^ — cxd. It is reassnring that (^, rather than (|^), is selected 
as the correct formnla, since dn in (H) goes to 2 for d —>• — cx). In addition (^2[) 
provides strong evidence for the existence of a nniqne fractal dimension dn = dh at 
all distances. As shown in the Appendix one cannot take for granted snch a relation, 
in particular for non-unitary theories coupled to quantum gravity, but based on (P2[ ) 
it is tempting to conjecture that, contrary to the situation for the so-called multi- 
critical branched polymers discussed in the Appendix, we will always have dn = dh 
in quantum gravity. 

We have emphasised the validity of (§) for c < 0, but have remained silent about 
the region 0 < c < 1. The reason is that for non-unitary conformal held theories we 
have operators with negative scaling dimensions. They will be dominant relative to 
the cosmological term, and it cannot be entirely ruled out that a formula like (|^) is 
correct for c > 0. Contrary to eq. (|^) it has a drastic dependence on c for c —> 1. 
This dependence has not been observed until now in the computer simulations which 
favour dn ~ 4. For c > 0 the numerical data are based on ordinary Monte Carlo 
simulations and it is not possible to perform the same high precision determination 
of dn as for c = —2. The best present data are not in good agreement with (|), but 
cannot falsify it either. 

Eqs. (|5^ and (|53D show a new and surprising scaling indicating that //r^ seems to 
be a universal, dimensionless variable. In the case c = 0 this is certainly reasonable 
since we know that the dimension of r is and it is natural to expect that 

a boundary has dimension This is not satisfied for {l)r due to the fractal 

structure of space-time, which implies that the boundary of a sphere of geodesic 
radius r is highly multi-connected with many microscopic loops. This manifests itself 
in the short distance cut-off needed in eq. ( 0 ) for n = 1. However, for moments 
n > 1, these microscopic loops play no role and we get ~ A^"-A (^g^e eq. 
(^)). If //r^ is a universal, dimensionless variable also for c < 0 where dn < 4, 
we reach the conclusion that the I appearing in the computer measurements cannot 
be identihed with the “Lionville” £, which in the continuum notation is believed to 
have the dressing dictated by: 

S( j - i), (93) 
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where a denotes the gravitational dressing associated with the cosmological term and 
the gravitational dressing of the boundary cosmological term is a/2. In the context 
of Liouville theory it has never been entirely clear how to derive this result from hrst 
principles since the actual boundary conditions of the matter helds for c < 0 have 
never been explicitly specified. In the numerical simulations the boundary is not 
fixed, but only characterised by being a (multi-connected) spherical shell of geodesic 
radius r. Probably free boundary conditions in Liouville theory come closest to the 
“experimental” set-up in the numerical simulations, and maybe the existence of 
scaling operators with negative dimensions in the non-unitary theories can spoil the 
relationship between the dimension of space-time and boundary shown in eq. (P3[ ) 
in the case of free boundary conditions, since these operators are dominant relative 
to the cosmological constant. 

We hope it will eventually be possible to understand this new, observed scaling 
of the boundary length from hrst principles. 
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Appendix 


The purpose of this appendix is to provide an example of a statistical ensemble to 
which one can assign both kind of fractal dimensions, dn and dh as discussed in 
sec. 2.3, but where they do not coincide. The model we have in mind is a model of 
so-called multi-critical branched polymers |^. It is dehned as a statistical ensemble 
of connected, planar graphs without any loops, i.e. connected planar tree-graphs. 
For such a tree-graph the weight is given by a fugacity factor e~^ for each link I and 
a branching factor f{nv) for each vertex v of order n^. BP denotes the class of 
planar tree graphs, the partition function is given by 


= E ^ n /(^-) n e 

G&BP y^G l£G 




( 94 ) 


where Sg is a symmetry factor for the graph G, such that a rooted branched polymer 
is counted only once. It is known that for a large class of positive functions / the 
hrst derivative of the partition function is given by 

= C 2 -ill -for (95) 

where ^2 and C 2 are non-universal constants. The geodesic distance between two 
vertices in a tree graph is dehned as the shortest link distance between the two 
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vertices. The two-point function, defined as in eq. (P^), except that two marked 
points are separated a geodesic distance r, can be calculated and is given by 




e 


k2 

5 


( 96 ) 


for /i —>■ where k 2 is again a non-universal constant. From (P3|) and (P^D we can 
hnd the partition function as well as the two-point function for constant volume, i.e. 
for a constant number of links, by a (discrete) Laplace transformation in /x: 


(r) ~ Ar-3/2 ^ rVTV ^ (^gg^ 

/o'\ 

The “spherical shell” of geodesic radius r for graphs of volume iV, n)^ (r), is dehned 
by analogue with (^) and counts the average number of vertices of distance r from 
an arbitrarily chosen vertex: 




G 


r] = 


( 2 ), 

N ^ 


NZ 


( 2 ) 

N 


~ r e 


-k2 r'^/N 


(99) 


Comparing (PTD-(P^) with (P^ we conclude that du = = 2 for this class of 

branched polymers. 

There exists a generalisation of the generic class of branched polymers considered 
here, which in many respects are related to the generic branched polymers as the 
multi-critical matrix models are related to pure two-dimensional gravity. They are 
called multi-critical branched polymers and as for multi-critical matrix models they 
are dehned by allowing certain negative weights for some of the orders of vertices. 
By hne-tuning of the weights one can obtain a new critical behaviour: 


Z^r^' = for /i 




( 100 ) 






( 101 ) 


where Cm, /im and km are non-universal constants and m > 2 (m = 2 corresponds 
to the generic branched polymer). From eq. (|101|) it follows that du = rn/{m — 1). 
However, inverse Laplace transformations lead to the following expressions (in the 
limit of large N) for the partition function with hxed volume N, 


zip'' ~ AI-2-'/™ X 


( 102 ) 


and the two-point function with hxed volume N, 

G^Y'\r) = n diaef^^G^^Yr) 

J—ioo 


' —200 
^200 


[ dfi (l-kmr{n-frmY + - 

J—ioo ^ ^ 


NZP(^ 


Co + Cl 


jYl-2/r; 


-h 0(r^)). 


(103) 
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From this equation we conclude that 


n^^\r) -^ for < r < ( 104 ) 

N rn 

This shows that from a formal point of view dh = 2, and one has the situation 
indicated in eq. (p^), with dh = 2 and dn = m/{m — 1). 
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N 

dh 

a 

C 


^min 

^max 

8192000 

3.5335(7) 

4.277(8) 

0.0429(1) 

1.2 

5 

105 

2048000 

3.5253(6) 

4.222(6) 

0.0443(1) 

1.1 

5 

60 

1024000 

3.5166(8) 

4.184(6) 

0.0456(2) 

1.5 

5 

45 

512000 

3.498(2) 

4.08(2) 

0.0489(4) 

1.0 

5 

45 

256000 

3.4932(6) 

4.061(4) 

0.0497(1) 

1.7 

5 

30 

128000 

3.483(3) 

4.020(8) 

0.0515(3) 

1.2 

5 

15 

64000 

3.464(1) 

3.945(6) 

0.0546(2) 

0.5 

5 

15 

8192000 

3.533(1) 

4.27(2) 

0.0430(2) 

1.2 

10 

105 

2048000 

3.525(1) 

4.21(1) 

0.0445(2) 

1.1 

10 

60 

1024000 

3.529(1) 

4.22(1) 

0.0452(3) 

1.1 

10 

50 

512000 

3.490(4) 

3.99(4) 

0.0507(8) 

0.9 

10 

45 

256000 

3.490(1) 

4.03(1) 

0.0504(3) 

1.4 

10 

30 

128000 

3.463(2) 

3.88(1) 

0.0555(3) 

3.6 

10 
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Table 1: The fractal dimension dh as determined from the small distance scaling of 
nisf{r) when r is the triangle distance. 
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N 

dh 

a 

C 


'^min 

^max 

512000 

3.593(5) 

0.480(9) 

1.72(2) 

1.0 

2 

10 

256000 

3.625(3) 

0.530(5) 

1.59(1) 

1.5 

2 

6 

128000 

3.623(3) 

0.530(4) 

1.59(1) 

1.7 

2 

5 

64000 

3.607(2) 

0.512(3) 

1.645(7) 

8.7 

2 

5 

512000 

3.556(6) 

0.40(1) 

1.91(3) 

1.2 

3 

12 

256000 

3.594(5) 

0.476(9) 

1.72(2) 

0.8 

3 

7 

128000 

3.584(4) 

0.465(8) 

1.75(2) 

1.0 

3 

6 

64000 

3.544(3) 

0.408(5) 

1.92(1) 

6.9 

3 

6 

512000 

3.525(7) 

0.31(2) 

2.10(5) 

1.0 

4 

13 

256000 

3.560(7) 

0.40(2) 

1.89(4) 

0.7 

4 

8 

128000 

3.535(6) 

0.36(1) 

2.00(3) 

1.0 

4 

7 


Table 2: The fractal dimension dh as determined from the small distance scaling of 
nisf{r) when r is the link distance. 


link dist. 

triangle dist. 

N 

dn 

a 

dn 

a 

3.602(20) 

0.50(15) 

3.560(16) 

4.60(35) 

256000 ~ 

128000 

3.610(14) 

0.50(7) 

3.552(10) 

4.50(20) 

128000 - 

64000 

3.621(10) 

0.54(5) 

3.538(6) 

4.30(10) 

64000 ~ 

32000 

3.634(4) 

0.55(1) 

3.520(8) 

4.20(10) 

32000 - 

16000 

3.608(12) 

0.50(7) 

3.555(5) 

4.55(15) 

256000 - 

64000 

3.612(7) 

0.50(4) 

3.544(4) 

4.40(10) 

128000 - 

32000 

3.630(8) 

0.55(3) 

3.532(8) 

4.30(15) 

64000 - 

16000 

3.610(8) 

0.50(5) 

3.549(7) 

4.45(15) 

256000 - 

32000 

3.618(5) 

0.52(5) 

3.538(12) 

4.35(20) 

128000 - 

16000 

3.575(8) 

0.30(5) 

3.573(8) 

5.0(2) 

256000 - 

16000* 


Table 3: The fractal dimension dfj as determined from collapsing n]^{r). Both 
dehnitions of distance are included. * refers to the graphical computation of du 
from the intersection points of the dfj{(i) curves shown in Fig. |] 
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link dist. 

triangle dist. 

N 

d}j Oj 

d}{ a 

3.571(2) 0.131(4) 

3.573(12) 0.138(30) 
3.574(3) 0.137(7) 

3.575(4) 0.141(7) 

3.578(7) 0.150(25) 

3.570(5) 4.93(7) 

3.574(20) 5.02(40) 
3.573(3) 4.97(6) 

3.576(3) 5.04(7) 

3.574(7) 5.00(17) 

256000 - 2000 

256000 - 4000 

256000 - 8000 

256000 - 16000 

256000 - 32000 


Table 4: The fractal dimension dn as determined from collapsing the intersection 
point of Ra,]\f{d) for c = —2. Both dehnitions of distance are included. 


n 

N 

2n 

a 

C 


^min 

^max 

2 

512000 

4.12(1) 

0.310(8) 

11.6(3) 

1.2 

1 

5 


256000 

4.117(7) 

0.312(4) 

11.6(2) 

4.1 

1 

5 


128000 

4.10(1) 

0.309(7) 

11.8(3) 

3.1 

1 

5 


64000 

4.122(8) 

0.317(4) 

11.4(2) 

4.4 

1 

4 


512000 

3.95(1) 

0.17(1) 

16.8(5) 

1.2 

2 

8 


256000 

3.98(1) 

0.21(1) 

15.5(4) 

2.8 

2 

6 


128000 

3.94(2) 

0.18(2) 

16.7(7) 

0.9 

2 

6 


64000 

3.93(1) 

0.17(1) 

17.2(5) 

2.6 

2 

5 

3 

512000 

6.09(2) 

0.326(9) 

0.65(3)xl0^ 

1.3 

1 

6 


256000 

6.14(1) 

0.345(5) 

0.59(2)x102 

2.5 

1 

5 


128000 

6.11(2) 

0.342(9) 

0.61(3)x102 

2.4 

1 

5 


64000 

6.14(2) 

0.348(6) 

0.58(2)x102 

2.5 

1 

4 


512000 

5.91(2) 

0.21(2) 

1.01(6)x102 

0.9 

2 

8 


256000 

5.94(2) 

0.23(1) 

0.93(5)x102 

2.8 

2 

6 


128000 

5.86(4) 

0.20(2) 

1.07(9)x102 

0.6 

2 

6 


64000 

5.83(3) 

0.19(1) 

1.12(6)x102 

2.5 

2 

5 

4 

512000 

8.02(2) 

0.325(9) 

0.50(3)xl0^ 

2.1 

1 

8 


256000 

8.18(2) 

0.373(7) 

0.36(2)x103 

1.4 

1 

5 


128000 

8.13(4) 

0.37(1) 

0.37(1)x103 

1.6 

1 

5 


64000 

8.01(2) 

0.336(6) 

0.47(2) xlO^ 

10 

1 

5 


512000 

7.81(3) 

0.22(2) 

0.85(7) xlO^ 

1.3 

2 

9 


256000 

7.90(4) 

0.26(2) 

0.68(6)x103 

2.5 

2 

6 


128000 

7.64(4) 

0.16(2) 

1.2(1) xlO^ 

1.9 

2 

7 


64000 

7.74(4) 

0.21(2) 

0.90(9)x103 

2.0 

2 

5 


Table 5: The exponent 2n as determined from the small distance scaling of {R)r,N- 


32 





n 

dfj Oj 

N 

2 

3.629 (33) 0.35 (20) 
3.616 (23) 0.20 (10) 
3.606 (18) 0.25 (10) 
3.588 (14) 0.15 (5) 

256000 - 64000 

128000 - 32000 

64000 - 16000 

32000 - 8000 

3 

3.654 (28) 0.40 (20) 
3.645 (25) 0.23 (17) 
3.636 (21) 0.27 (10) 
3.621 (20) 0.20 (10) 

256000 - 64000 

128000 - 32000 

64000 - 16000 

32000 - 8000 

4 

3.662 (40) 0.40 (30) 
3.648 (32) 0.20 (20) 
3.641 (33) 0.25 (15) 
3.626 (30) 0.20 (15) 

256000 - 64000 

128000 - 32000 

64000 - 16000 

32000 - 8000 


Table 6: The fractal dimension du as determined from collapsing {l^) 
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Figure 4: (a) w{n,k) for n = 16000 plotted from Eq. (^Tj) and by application of 
our program routine 10® times. (6) The deviation Aw{n,k) = (ta(n, A:)theoretical — 
w{n, Ai)measured)/w(n, Ai)theoretical usiug the Same parameters. 
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1.6 


16000 
32000 
64000 
. 128000 
256000 


-2 


-1 
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Figure 5: (a) The small x behaviour of the logarithmic derivative of n]^{r) where r 
is the triangle distance. We use dn = 3.50 and a = 4.11. (b) The same as in (a) bnt 
now njqij-) is dehned in terms of link distance, dn = 3.62 and a = 0.52. 
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Figure 6: (a) The ni^{r) distributions defined in terms of triangle distance rescaled 
according to Eq. (|30|) using dn = 3.56 and a = 4.50. (6) The same as in (a) but 
now by setting a = 0. 
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Figure 7: (a) The n]sf{r) distributions defined in terms of link distance rescaled 
according to Eq. (|30|) using dn = 3.60 and a = 0.50. (6) The same as in (a) but 
now by setting a = 0. 
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3.75 
3.7 
3.65 

3.6 
3.55 

3.5 
3.45 

4 5 6 

a 

3.75 

3.7 
3.65 

-c 3.6 
3.55 
3.5 
3.45 

Figure 8: (a) dff(a) computed from collapsing the njq{r) distributions defined in 
terms of triangle distance. We collapse the distributions pairwise and in the plot we 
indicate the largest lattice. For clarity, errors computed from are displayed only 
for the largest lattice. (6) The same as in (a) but now nN{r) is dehned in terms of 
link distance. 
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dH(a) 


Figure 9: (a) dn computed from Xmm by collapsing the njq{r) distributions defined 
in terms of triangle distance. Each point on the graph is dnia) (see Fig. |]) and is 
plotted versus the corresponding xLini^)- (^) The same as in (a) but now nN{r) is 
dehned in terms of link distance. 
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Figure 10: (a) The functions -Ra,N(d), defined in terms of link distance, for N = 

2K, 4:K ,..., 256K and a = 0.130. (6) Same as in (a) but when Ra^N^d) is defined in 
terms of triangle distance and a = 5.0. 
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dh 



dh 


Figure 11: (a) defined by Eq. (|46D for two sets of iVj’s when Ra^N{d) is 

defined in terms of link distance. (6) Same as in (a) but when Ra^N^d) is defined in 
terms of triangle distance. 
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Figure 12: The small x behaviour of the logarithmic derivative for {l"')r,N for n = 2, 
3, 4, iV = 16000-512000 using a = 0.32 and dn = 3.63, 3.65 and 3.66 respectively. 
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Figure 13: (a) The {P)r,N distributions defined in terms of link distance rescaled 
according to Eq. (^) using du = 3.63 and a = 0.35. (6) The same as in (a) for 
{l^)r,N, dn = 3.65 and a = 0.4. (c) The same as in (a) for dn = 3.66 and 

a = 0.4. 
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Figure 15: The consequences of scaling using the assumption dim[iV] = dim[/^]: (a) 
The rescaled according to Eq. (^) {f‘)r,N- {b) The loop length distribution function 
using y = dn = 3.58. 
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Figure 16: (a) The loop length distribution p^{r,l) for N = 64000,128000 and 
256000. (fc) The same as in (a) for N = 64000 compared to the corresponding 
observable for pure gravity (c = 0). The dashed line is a £t to Eq. (^). The data 
for c = 0 was taken from [^. (c) The same as in (a), but for (approximately) hxed 
and du = 3.58. 
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